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The  goal  of  this  MURI  was  to  detect  magnetic  resonance  from  an  individual  proton  using  magnetic  resonance  force  microscopy  (MRFM). 
Pushing  MRFM  to  single-proton  sensitivity  requires  meeting  three  experimental  challenges:  (1)  Fabricating  nano  magnets  which  produce  a 
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Scientific  Progress 

The  goal  of  this  MURI  was  to  detect  magnetic  resonance  from  an  individual  proton  using  magnetic  resonance  force  microscopy 
(MRFM).  Pushing  MRFM  to  single-proton  sensitivity  requires  meeting  three  experimental  challenges: 

1 .  Fabricating  nanomagnets  which  produce  a  large  magnetic  field  gradient  and  yet  have  only  a  thin  damage  layer.  An  MRFM 
experiment  may  be  performed  in  either  a  "sampleon-  cantilever”  configuration,  with  the  sample  affixed  to  the  cantilever  and 
brought  near  an  independently  fabricated  nanomagnet,  or  in  a  ”magnet-on-cantilever”  configuration,  with  the  nanomagnet 
affixed  to  the  cantilever  and  brought  near  a  thinfilm  sample.  To  allow  for  imaging  the  widest  range  of  samples  —  including 
semiconductor  devices  and  flash-frozen  biological  samples  —  we  focused  attention  on  preparing  nanomagnets  affixed  to  the 
leading  edges  of  attonewton-sensitivity  cantilevers  [1;  2], 

2.  Understanding  and  mitigating  the  large  non-contact  friction  (and  frequency  noise)  typically  experienced  by  an 
attonewton-sensitivity  cantilever  near  a  surface.  In  the  highestsensitivity  magnetic  resonance  force  microscope  experiments 
carried  out  to  date,  the  cantilever  force  sensitivity  has  been  limited  not  by  intrinsic  dissipation  in  the  cantilever,  but  by  spurious 
force  (and  frequency)  fluctuations  experienced  by  the  cantilever  near  a  surface.  We  have  developed  a  first-principles 
understanding  of  this  surface  noise  [3-8]  in  order  to  minimize  it  through  improved  cantilever  design  [1 ;  2]  and  sample 
preparation. 

3.  Devising  a  suitable  spin  detection  protocol.  The  signals  in  an  MRFM  experiment  are  small.  A  central  challenge  in  MRFM  is 
devising  spin-detection  protocols  which  (1)  perform  well  in  the  presence  of  surface  noise  and  (2)  are  immune  from  artifacts 
caused  by  the  microwave  or  radiowave  perturbations  applied  to  manipulate  spin  magnetization.  The  “i-OSCAR”  spin-detection 
protocol  invented  by  Dan  Rugar  and  co-workers  to  detect  single-electron  magnetic  resonance  via  MRFM  met  these  exacting 
demands  [9],  but  required  samples  with  extraordinarily-long  spin  relaxation  times  to  achieve  single  spin  sensitivity.  To  extend 
single-electron  spin  detection  to  scientifically  and  technologically  interesting  samples  clearly  requires  new  spin-detection 
protocols  [10;  11].  For  single-proton  spin  detection,  we  have  tried  to  determine  whether  force  detection  or  force-gradient 
detection  is  more  likely  to  succeed  in  generating  a  distinguishable  signal. 

Pushing  MRFM  to  single-proton  sensitivity  also  requires  meeting  a  number  of  theoretical  challenges:  Devising 
image-reconstruction  algorithms  that  work  in  the  single-spin  limit,  when  the  signal-to-  noise  ratio  is  low  and  when  the 
experimental  knowledge  of  the  shape  and  magnetization  of  the  magnetic  tip  is  uncertain.  2.  Inventing  an  approach  for 
simulating  large  numbers  of  spins,  so  that  we  may  obtain  a  first-principles  understanding  of  spin  diffusion  and  dynamic  nuclear 
polarization  in  the  presence  of  the  large  magnetic  field  gradients  present  in  an  MRFM  experiment. 
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of  Magnet-Tipped  Attoneivton-Sensitivity  Cantilevers,  American  Vacuum  Society 
58th  International  Symposium  and  Exhibition;  Nashville,  Tennessee;  October 
30  -  November  4,  2011,  [url], 

3.  L.  Harrell,  E.  Moore,  S.  Lee,  S.  Hickman,  and  J.  Marohn,  Parametric  Amplifica¬ 
tion  Protocol  for  Frequency-Modidated  Magnetic  Resonance  Force  Microscopy  Sig¬ 
nals,  American  Physical  Society  March  Meeting;  Dallas,  Texas;  March  21  -  25, 
2011,  [url], 

4.  J.  A.  Marohn,  E.  Moore,  and  J.  Longenecker,  Rapid  Serial  Prototyping  and  Analysis 
of  Nanomagnet-Tipped  Attonewton-Sensitivity  Cantilevers  for  Magnetic  Resonance 
Force  Microscopy,  American  Physical  Society  March  Meeting;  Dallas,  Texas; 
March  21  -  25,  2011,  [url]. 
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5.  E.  W.  Moore,  S.  Lee,  S.  A.  Hickman,  J.  G.  Longenecker,  and  J.  A.  Marohn, 
Cantilever  torque  magnetometry  studies  of  the  in-plane  to  out-of-plane  transition 
in  a  single  nickel  magnetic  nanorod,  American  Physical  Society  March  Meeting; 
Dallas,  Texas;  March  21  -  25,  2011,  [url], 

6.  J.  G.  Longenecker,  S.  A.  Hickman,  E.  W.  Moore,  S.-G.  Lee,  S.  J.  Wright,  L.  E.  Har¬ 
rell,  and  J.  A.  Marohn,  Fabrication  of  Integrated  Nanomagnets  Overhanging  Batch- 
Fabricated  Attonewton-Sensitivity  Cantilevers,  57th  International  Symposium  and 
Exhibition;  Albuquerque,  New  Mexico;  October  17  -  22,  2010,  [url], 

7.  E.  W.  Moore,  S.-G.  Lee,  S.  A.  Hickman,  S.  J.  W.  andLee  E.  Harrell,  J.  G.  Longe¬ 
necker,  P.  P.  Borbat,  J.  H.  Freed,  and  J.  A.  Marohn,  High  sensitivity  electron  spin 
magnetic  resonance  force  microscopy  for  labeled  biological  samples,  American  Phys¬ 
ical  Society  March  Meeting;  Portland,  Oregon;  March  15  -  19,  2010,  [url], 

8.  S.-G.  Lee,  E.  W.  Moore,  S.  A.  Hickman,  L.  E.  Harrell,  and  J.  A.  Marohn, 
Non-degenerate  parametric  amplification  used  for  surface  noise  evasion  in  scanned 
probe  microscopy,  American  Physical  Society  March  Meeting;  Portland,  Oregon; 
March  15  -  19,  2010,  [url]. 

9.  J.  G.  Longenecker,  S.  Hickman,  L.  Harrell,  and  J.  Marohn,  Fabrication  of  over¬ 
hanging  magnet-tipped  cantilevers  for  nanoscale  scanned-probe  magnetic  resonance, 
American  Vacuum  Society  56th  International  Symposium  and  Exhibition;  San 
Jose,  California;  November  8  -  13,  2009,  [url]. 

10.  S.-G.  Lee,  New  Physics  for  Coupling  Spin  Magnetization  to  a  Mechanical  Oscil¬ 
lator,  Molecular  Imaging  2009:  Routes  to  Three-Dimensional  Imaging  of  Single 
Molecules;  Ithaca,  New  York;  August  9  -  13,  2009,  [url], 

11.  S.  A.  Hickman,  Batch  Fabrication  of  Overhanging  sub-100  nm  Diameter  Nickel 
Nanorods  on  Attonewton  Sensitivity  Silicon  Cantilevers,  Molecular  Imaging  2009: 
Routes  to  Three-Dimensional  Imaging  of  Single  Molecules;  Ithaca,  New  York; 
August  9  -  13,  2009,  [url], 

2.3  Intellectual  property 

•  Number  of  Patents  Disclosed:  0 

•  Number  of  Patents  Awarded:  0 

2.4  Technology  transfer 

Any  specific  interactions  or  developments  which  zvoidd  constitute  technology  transfer  of  the 
research  residts.  Examples  include  interaction  with  other  DOD  scientists,  interactions  with 
industry,  initiation  of  a  start-  up  company  based  on  research  results  or  transfer  of  information 
which  might  impact  the  development  of  products. 
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•  SC  Solutions.  The  Cornell  team,  in  collaboration  with  Dr.  Doran  Smith 
of  the  U.S.  Army  Research  Laboratory,  previously  worked  together  in  an 
ARO-funded  Phase  II  Small  Business  Technology  Transfer  Research  Grant 
with  SC  Solutions  (Sunnyvale,  CA)  and  Cryolndustries  of  America  (Manch¬ 
ester,  NH).  This  collaboration  previously  resulted  in  two  joint  publications 
[52;  53].  During  the  reporting  period,  SC  Solutions  delivered  an  FPGA-based 
cantilever  controller  to  the  marketplace.  Product  details  can  be  found  at 
http:/ / www.scsolutions.com/ public/ product/ mrfm-controllers.html. 

•  IBM.  The  Hero-team's  interactions  with  Dan  Rugar  and  coworkers  led  to  acknowl¬ 
edgements  of  these  contributions  in  his  recent  paper  [54],  Graduate  student  Jonilyn 
Longenecker  and  postdoc  Lei  Chen  of  Cornell  University  spent  a  week  at  IBM 
Almaden  Research  Center  to  execute  a  joint  research  project  with  Dr.  John  Mamin 
and  Dr.  Dan  Rugar  (Almaden,  California;  November  10  -  18,  2011).  They  had  previ¬ 
ously  visited,  along  with  Marohn,  the  laboratory  of  Daniel  Rugar  at  IBM  Almaden 
to  learn  advanced  experimental  techniques  in  magnetic  resonance  force  microscopy 
(Almaden,  California;  October  12  -  13,  2009).  The  Cornell/IBM  collaboration  will 
result  in  one  joint  publication  [6]. 

•  USMA.  Marohn's  team  at  Cornell  hosted  and  carried  out  joint  research  with  Prof. 
Lee  Harrell  of  the  U.S.  Military  Academy  at  West  Point  (June  2010  and  June  2011). 
The  Cornell/USMA  collaboration  resulted  in  three  joint  publications  [12-14], 

•  USARL.  Visited,  along  with  two  graduate  students,  the  laboratory  of  Dr.  Doran 
Smith  at  the  U.S.  Army  Research  Laboratory  to  learn  advanced  experimental  tech¬ 
niques  in  magnetic  resonance  force  microscopy  (Adelphi,  Maryland;  April  24, 2010). 
Supplied  attonewton-sensitivity  cantilevers  to  Smith's  team  and  interacted  with 
him  via  email  biweekly  (approximately)  during  the  reporting  period.  The  Cor- 
nell/USARL  collaboration  resulted  in  one  joint  publication  [8]. 

•  NIH.  Visited,  along  with  two  graduate  students,  the  laboratory  of  Dr.  Sriram  Subra- 
maniam  at  the  National  Institutes  of  Health  to  learn  about  cryo-electron  microscopy 
and  preparing  biological  samples  for  study  in  vacuum  at  cryogenic  temperatures 
(Bethesda,  Maryland;  April  21,  2010). 

•  Interacted  with  Dr.  Doran  Smith  (US  ARL)  and  Drs.  John  Mamin  and  Dan  Rugar 
(IBM  Almaden)  during  Army  Research  Office  Multi-University  Research  Institute 
(ARO/MURI)  program  review  meetings. 

-  Almaden,  California;  September  29  -  30,  2011. 

-  Washington,  DC;  April  22  -  23,  2010. 
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3  Scientific  progress  and  accomplishments 

3.1  Introduction 

The  goal  of  this  MURI  was  to  detect  magnetic  resonance  from  an  individual  proton  using 
magnetic  resonance  force  microscopy  (MRFM).  Pushing  MRFM  to  single-proton  sensi¬ 
tivity  requires  meeting  three  experimental  challenges: 

1.  Fabricating  nanomagnets  which  produce  a  large  magnetic  field  gradient  and  yet  have  only 
a  thin  damage  layer.  An  MRFM  experiment  may  be  performed  in  either  a  "sample- 
on-cantilever"  configuration,  with  the  sample  affixed  to  the  cantilever  and  brought 
near  an  independently  fabricated  nanomagnet,  or  in  a  "magnet-on-cantilever"  con¬ 
figuration,  with  the  nanomagnet  affixed  to  the  cantilever  and  brought  near  a  thin- 
film  sample.  To  allow  for  imaging  the  widest  range  of  samples  —  including  semi¬ 
conductor  devices  and  flash-frozen  biological  samples  —  we  focused  attention  on 
preparing  nanomagnets  affixed  to  the  leading  edges  of  attonewton-sensitivity  can¬ 
tilevers  [1;  2], 

2.  Understanding  and  mitigating  the  large  non-contact  friction  (and  frequency  noise)  typi¬ 
cally  experienced  by  an  attonewton-sensitivity  cantilever  near  a  surface.  In  the  highest- 
sensitivity  magnetic  resonance  force  microscope  experiments  carried  out  to  date, 
the  cantilever  force  sensitivity  has  been  limited  not  by  intrinsic  dissipation  in  the 
cantilever,  but  by  spurious  force  (and  frequency)  fluctuations  experienced  by  the 
cantilever  near  a  surface.  We  have  developed  a  first-principles  understanding  of 
this  surface  noise  [3-8]  in  order  to  minimize  it  through  improved  cantilever  design 
[1;  2]  and  sample  preparation. 

3.  Devising  a  suitable  spin  detection  protocol.  The  signals  in  an  MRFM  experiment  are 
small.  A  central  challenge  in  MRFM  is  devising  spin-detection  protocols  which 
(1)  perform  well  in  the  presence  of  surface  noise  and  (2)  are  immune  from  arti¬ 
facts  caused  by  the  microwave  or  radiowave  perturbations  applied  to  manipulate 
spin  magnetization.  The  "i-OSCAR"  spin-detection  protocol  invented  by  Dan  Rugar 
and  co-workers  to  detect  single-electron  magnetic  resonance  via  MRFM  met  these 
exacting  demands  [9],  but  required  samples  with  extraordinarily-long  spin  relax¬ 
ation  times  to  achieve  single  spin  sensitivity.  To  extend  single-electron  spin  detec¬ 
tion  to  scientifically  and  technologically  interesting  samples  clearly  requires  new 
spin-detection  protocols  [10;  11].  For  single-proton  spin  detection,  we  have  tried 
to  determine  whether  force  detection  or  force-gradient  detection  is  more  likely  to 
succeed  in  generating  a  distinguishable  signal. 

Pushing  MRFM  to  single-proton  sensitivity  also  requires  meeting  a  number  of  theo¬ 
retical  challenges: 

1.  Devising  image-reconstruction  algorithms  that  work  in  the  single-spin  limit,  when  the 
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signal-to-  noise  ratio  is  low  and  when  the  experimental  knowledge  of  the  shape  and 
magnetization  of  the  magnetic  tip  is  uncertain. 

2.  Inventing  an  approach  for  simulating  large  numbers  of  spins,  so  that  we  may  obtain  a 
first-principles  understanding  of  spin  diffusion  and  dynamic  nuclear  polarization  in 
the  presence  of  the  large  magnetic  field  gradients  present  in  an  MRFM  experiment. 
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3.2  Accomplishments 

This  MURI's  funding  began  on  2005/03/01,  ended  on  2010/09/30,  and  was  continued 
under  a  no-cost  extension  until  2011/09/30.  Here  we  summarize  progress  made  during 
the  entire  project,  highlighting  progress  made  during  the  final  reporting  period  from 
2009/08/01  to  2011/09/30. 

3.2.1  Fabricating  nanomagnets  on  attonewton-sensitivity  cantilevers 

Addressing  this  challenge  was  the  Marohn  team's  main  contribution  to  the  MURI  effort. 

Spin  sensitivity  in  a  magnetic  resonance  force  microscope  measurement  is  propor¬ 
tional  to  the  cantilever's  force  sensitivity  and  inversely  proportional  to  the  available  mag¬ 
netic  field  gradient.  Achieving  a  large  magnetic  field  gradient  requires  fabricating  mag¬ 
nets  with  a  small  diameter  and  bringing  the  sample  as  close  as  possible  to  the  magnet. 
It  is  thus  essential  that  the  cantilever  maintain  its  force  sensitivity  at  small  cantilever- 
sample  separations  (in  magnet-on-cantilever  experiments)  or  at  small  cantilever-magnet 
separations  (in  sample-on-cantilever  experiments). 


3.2.1. 1  Radio  frequency  cantilevers 

In  the  first  half  of  the  MURI,  the  Marohn  team  developed  both  audiofrequency 
and  radiofrequency  cantilevers  for  MRFM.  Audio  frequency  cantilevers  had  exhibited 
attonewton  sensitivity  [12;  13]  provided  they  were  made  sufficiently  thin,  narrow,  and 
long,  and  were  cooled  to  temperatures  of  4.2  kelvin  or  below.  Radio  frequency  cantilevers 
had  also  been  developed  [14;  15],  but  at  the  time  the  MURI  began  the  force  sensitivity  of 
rf  cantilevers  was  inferior  to  the  force  sensitivity  of  audio  frequency  cantilevers  [16]. 

Reasoning  as  follows,  we  nevertheless  expected  rf  cantilevers  to  better  maintain  their 
force  sensitivity  at  small  tip-sample  separations.  According  to  the  fluctuation-dissipation 
theorem,  the  force  sensitivity  of  a  cantilever  is  proportional  to  temperature  and  the  dis¬ 
sipation  (e.g.,  damping)  experienced  by  the  cantilever.  A  major  source  of  cantilever 
damping  near  a  surface  arises  from  charges  on  the  cantilever  interacting  with  surface 
electric  field  fluctuations  in  the  sample  [17],  with  the  cantilever  dissipation  proportional 
to  the  spectral  density  of  electric  field  fluctuations  at  the  cantilever  resonance  frequency. 
Over  polymeric  samples  these  fluctuating  electric  fields  arise  from  thermal  dielectric  fluc¬ 
tuations  [3;  4],  The  dielectric  response  of  most  polymers,  including  biopolymers,  falls  off 
at  frequencies  above  a  few  hundred  kHz.  One  would  therefore  expect  a  cantilever  with 
a  resonance  frequency  in  the  MHz  range  to  exhibit  much  less  dissipation  near  a  polymer 
surface  than  a  cantilever  with  a  resonance  frequency  in  the  kHz  range. 

When  we  began  the  MURI,  there  were  many  examples  of  radiofrequency  cantilevers, 
but  the  great  majority  were  doubly-clamped  beams  and  therefore  unsuitable  for  scanned 
probe  experiments.  A  second  challenge,  when  we  began,  was  devising  a  scheme  for 
detecting  the  motion  of  the  rf  cantilever.  Laser  interferometry  has  been  used,  but  did  not 
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have  the  sensitivity  to  detect  the  thermal  Brownian  motion  in  the  (high  spring  constant) 
rf  cantilever.  The  most  commonly  used  detection  scheme  was  magnetomotive  detection. 
This  method  required  the  cantilever  beam  to  be  attached  at  each  end  to  a  substrate  and 
was  therefore  unsuitable  for  scanned  probe  experiments. 

In  previous  reports  we  described  in  detail  our  (unpublished)  work  developing  silicon 
rf  cantilevers  by  top-down  nanofabrication  for  MRFM  experiments,  which  we  briefly 
summarize  here.  We  began  by  fabricating  doubly  clamped  beams  and  observing  their 
motion  via  magneto-motive  detection,  demonstrating  that  operating  at  higher  applied 
magnetic  field  for  increased  sensitivity  also  led  to  more  readout-induced  damping  due  to 
current  fluctuations  in  the  beam  interacting  with  the  applied  magnetic  field.  In  parallel, 
we  invented  and  developed  a  method  for  fabricating  an  rf  cantilever  with  an  integrated 
tunnel  sensor  for  reading  out  cantilever  motion;  we  succeeded  in  making  tunnel  sensors 
by  opening  up  a  gap  in  a  metal-film  constriction  using  electromigration,  although  not  on 
a  tunnel  sensor  integrated  with  a  cantilever.  Finally,  we  succeeded  in  fabricating  singly- 
clamped  cantilevers  with  integrated  (but  untested)  capacitance  sensors;  these  cantilevers 
were  free-standing  and  in  principle  suitable  for  scanned  probe  experiments. 

At  the  end  of  Phase  I  of  our  MURI,  we  made  the  difficult  decision  to  abandon  our 
team's  rf  cantilever  work  due  to  budget  constraints. 


3.2.1.2  Audio  frequency  cantilevers 

When  we  began  the  MURI,  both  ESR  and  NMR  had  been  detected  in  a  magnet-on- 
cantilever  MRFM  experiment.  The  magnetic  tips  in  these  experiments  were  prepared 
by  affixing  a  few-micron  diameter  magnetic  particle  to  the  cantilever's  leading  edge  by 
hand,  with  further  focused  ion  beam  milling  used  to  prepare  submicron  tips. 

Our  team  had  previously  demonstrated  patterning  of  100  nm  scale  nickel  magnets  on 
attonewton  sensitivity  cantilevers  [18],  but  the  yield  was  low  and  alignment  challenges 
left  us  unable  at  that  time  to  prepare  tips  close  enough  to  the  cantilever's  leading  edge  for 
an  MRFM  experiment.  Sean  Garner  in  Marohn's  group  invented  and  demonstrated  the 
idea  of  using  anisotropic  wet  etching  of  silicon  to  prepare  magnets  overhanging  a  silicon 
edge.  Garner  demonstrated  the  process  with  nickel  magnets  defined  using  optical  lithog¬ 
raphy  and  this  work  was  extended  by  Marohn  to  magnets  defined  by  e-beam  lithography. 
Unfortunately,  the  process  required  silicon-on-insulator  wafers  with  a  (111) -oriented  sil¬ 
icon  device  layer;  these  became  impossible  to  obtain  in  2005,  and  the  process  could  not  be 
used  to  make  magnetic-tipped  cantilevers. 

With  MURI  funding,  Hickman  at  Cornell  was  able  to  attack  the  challenge  of  pro¬ 
ducing  nanomagnetic  tips  overhanging  the  leading  edge  of  a  cantilever.  This  work  was 
summarized  in  detail  in  previous  reports.  The  key  innovation  was  identifying  and  opti¬ 
mizing  a  dry  isotropic  silicon  etch  that  did  not  etch  nickel  or  cobalt.  It  should  be  appre¬ 
ciated  that  harnessing  this  innovation  to  produce  MRFM  cantilevers  required  a  tremen¬ 
dous  amount  of  process  engineering  and  integration  —  the  final  cantilever  fabrication 
sequence  required  carrying  out  over  40  processing  steps,  mastering  18  nanofabrication 
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instruments,  and  required  approximately  three  weeks  of  full-time  effort  in  the  Cornell 
NanoScale  Science  and  Technology  Facility  This  reporting  period  we  published  the 
results  of  Hickman's  work  in  his  PhD  thesis  [19]  and  in  a  long  paper  in  ACS  Nano  [1].  In 
Ref.  1  we  demonstrated  using  one  of  Hickman's  nickel-tipped  cantilevers  to  detect  elec¬ 
tron  spin  resonance.  Moreover,  these  cantilevers  exhibited  a  record  low  surface-induced 
dissipation  (for  an  audio  frequency  cantilever)  at  small  tip-sample  separations. 

These  findings  established  that  (1)  e-beam  defined  nanomagnet  tips  could  be  fabri¬ 
cated  en  batch  overhanging  the  leading  edges  of  attonewton-sensitivity  cantilevers  and 
(2)  these  nanomagnet-tipped  cantilevers  could  be  used  to  detect  MRFM.  The  yield  was 
extremely  low,  however.  Furthermore,  the  agreement  between  the  observed  and  the 
predicted  field-dependence  of  the  ESR-MRFM  signal  was  poor,  indicating  that  the  tip 
was  incompletely  and  non-uniformly  magnetized.  A  subsequent  student  on  Marohn's 
team,  Jonilyn  Longenecker,  investigated  the  cause  of  the  low  yield.  The  primary  problem 
was  that  the  BOSCH  etch  used  to  remove  the  handle-wafer  silicon  beneath  the  cantilever 
involved  high  temperatures  which  promoted  chemical  reactions  between  the  nickel  or 
cobalt  of  the  magnetic  tip  and  the  silicon  in  the  underlying  wafer.  This  finding,  and  our 
team's  studies  of  barrier  layers,  were  likewise  detailed  in  prior  reports. 

During  this  reporting  period,  Longenecker  developed  a  method  for  fabricating  e-beam 
defined  nickel  and  cobalt  tips  overhanging  the  leading  edges  of  attonewton-sensitivity 
cantilevers.  Her  insight  was  to  decouple  the  magnet  fabrication  from  the  cantilever  fab¬ 
rication  to  avoid  exposing  the  magnets  to  high-temperature  processing.  To  achieve  this 
decoupling,  magnets  are  prepared  via  e-beam  lithography,  evaporation,  and  liftoff  on  sep¬ 
arate  micron-scale  "chips".  These  magnet-tipped  chips  are  released,  one  at  a  time,  from 
their  supporting  substrate  using  focused  ion  beam  milling  and  attached  to  a  cantilever  via 
focused  ion  beam  deposition.  Although  the  chips  are  affixed  to  the  cantilevers  serially,  the 
attachment  process  takes  only  a  few  hours  and  has  excellent  yield.  Both  the  magnets  and 
the  cantilevers  can  be  prepared  beforehand  via  batch  fabrication  in  high  yield.  A  paper 
describing  this  approach  to  fabricating  nickel  nanomagnet  tips  was  published  during  this 
reporting  period  in  Journal  of  Vacuum  Science  and  Technology  B  [20]. 

This  semi-serial  approach  to  producing  magnet-tipped  cantilevers  greatly  facilitates 
the  optimization  of  magnet  dimensions  and  the  analysis  of  identically-processed  mag¬ 
netic  material  by  X-ray  photoelectron  spectroscopy  (XPS)  and  SQUID  magnetometry.  It 
also  facilitates  the  analysis  of  individual  magnets  by  high-resolution  electron  microscopy 
and  electron  spectroscopy.  The  ability  to  analyze  nanomagnets  and  magnetic  films  in 
chemical  detail  at  nanometer  spatial  resolution  is  an  essential  step  in  producing  magnets 
for  single-proton  experiments,  where  less  than  5  nm  of  damage  at  the  magnet's  leading 
edge  is  required. 

In  Ref.  20  the  magnetic  properties  of  an  individual  cobalt  nanomagnet  on  a  can¬ 
tilever  were  assessed  using  frequency-shift  cantilever  magnetometry  [21-23].  To  suc¬ 
cessfully  fabricate  cobalt  nanomagnets,  the  process  of  Ref.  20  was  modified  to  elim¬ 
inate  high-temperature  processing  steps.  In  2011/11,  these  tips  were  used  to  detect 
proton  magnetic  resonance  in  collaboration  with  Dan  Rugar  and  John  Mamin  of  the  IBM 
Almaden  Research  Center.  Tips  were  prepared  at  Cornell,  flown  to  San  Jose  by  Lone- 
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Figure  1:  Cornell/IBM  magnetic  resonance  force  microscope  NMR  experiment,  (a)  Schematic  of 
the  experiment,  (b)  The  magnet-tipped  cantilever  used  in  the  experiment.  1.  Top-view  image  of 
the  magnet-tipped  chip  attached  to  an  underlying  cantilever  by  ion-assisted  platinum  deposition; 
three  rectangular  platinum  adhesion  patches  can  be  seen  on  the  top  side  of  the  chip.  The  over¬ 
hanging  cobalt  magnet  is  attached  to  a  3  pm  long  finger  at  the  leading  edge  of  the  chip.  Scale  bar 
=  5  pm.  2.  Angled  image  of  the  overhanging  cobalt  nanomagnet,  acquired  before  it  was  attached 
to  a  cantilever.  The  magnet  is  225  ±  15  nm  wide,  1494  ±  15  nm  long,  and  79  ±  4  nm  thick.  There 
is  a  4  nm  titanium  layer  under  the  magnet  to  promote  adhesion  to  the  silicon  substrate,  as  well 
as  an  8  nm  platinum  capping  layer  to  mitigate  oxidation.  Scale  bar  =  200  nm.  3.  Top-view  image 
of  a  custom-fabricated  200  pm  long  cantilever,  drawn  from  the  same  batch  as  the  cantilever  used 
in  this  experiment.  Scale  bar  =  20  pm.  (c)  Magnetic  resonance  signal  of  protons  in  a  40  nm  thick 
polystyrene  film.  Figure  and  caption  adapted  from  Ref.  20. 


necker  and  Marohn-group  postdoc  Lei  Chen,  and  inserted  into  the  IBM  magnetic  res¬ 
onance  force  microscope.  Proton  magnetic  resonance  experiments  were  carried  out  by 
Mamin  and  Rugar  in  collaboration  with  the  Cornell  team.  The  cobalt-tipped  cantilevers 
were  used  to  detect  stochastic  proton  magnetization  from  a  polystyrene  film  spun-cast 
over  a  microwire. 

See  Fig.  l(a-b).  The  experiments  employed  a  cobalt  magnet  with  cross  section  225  nm 
by  79nm;  note  how  the  magnet  extends  past  the  leading  edge  of  the  attonewton-sensitivity 
cantilever.  The  cantilever  was  centered  over  a  1  pm-wide  microwire  coated  with  40  nm 
of  polystyrene.  An  external  2.63  T  magnetic  field  was  applied  in  the  direction  of  the 
long  axis  of  the  cantilever  as  shown.  A  frequency-chirped  rf  waveform  (peak-to-peak 
deviation  2  MFIz)  was  applied  to  the  microwire  to  cyclically  invert  the  sample's  proton 
magnetization  at  twice  the  mechanical  resonance  frequency  fc  of  the  cantilever.  The  can- 
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tilever  motion  was  measured  with  a  laser  interferometer.  The  interferometer  output  was 
sent  to  a  lock-in  amplifier  that  inferred  the  Fourier  components  of  the  cantilever  oscil¬ 
lation  moving  in-phase  and  out-of-phase  with  the  rf  modulation.  Because  the  sample's 
proton  magnetization  had  a  coherence  time  of  ~  0.1  s,  the  sample's  magnetization  inter¬ 
acting  with  the  tip's  magnetic  field  gradient  led  to  spin-induced  force  fluctuations  acting 
on  the  cantilever.  A  spin  variance  signal  was  obtained  by  comparing  the  variance  of  the 
cantilever  displacement  observed  in  the  in-phase  and  out-of-phase  lock-in  channels. 

The  resulting  spin  variance  signal,  cvpin,  is  plotted  in  Fig.  1(c)  as  a  function  of  rf  fre¬ 
quency  at  three  tip-sample  separations.  The  peak  spin  signal  arises  from  a  root-mean- 
square  force  of  only  55  aN.  The  acquisition  time  was  12  to  17  minutes  per  data  point.  At 
this  acquisition  time  the  minimum  detectable  magnetic  moment  corresponds  to  approxi¬ 
mately  500  fully  polarized  protons.  The  magnetic  integrity  of  the  tip  at  the  nanoscale  was 
assessed  by  comparing  the  spin  signal  to  simulations  carried  out  using  different  models  of 
tip  damage  discussed  below.  The  tip-field  gradient  calculated  for  the  various  tip  models 
was  dBfP  jdz  =  4.5  to  5.4  MT  m-1. 

The  model  which  best  described  the  data  approximated  the  tip  as  fully  magnetized 
with  a  51  nm  thick  magnetically  inactive  layer  at  the  leading  edge.  The  simulated  signal 
is  shown  as  solid  blue  lines  in  Fig.  1(c);  the  agreement  between  experiment  and  simulation 
is  excellent.  The  magnetic  field  gradient  produced  by  our  cobalt  tip  is  7.6  to  9.1-fold  larger 
than  the  best  tip  gradient  achieved  in  any  previous  magnet-on-cantilever  MRFM  experi¬ 
ment  [24]  —  an  extremely  exciting  advance.  Moreover,  our  tip's  gradient  is  comparable  to 
the  4.2  MT  m-1  field  gradient  produced  by  the  Fe70Co3o  pillar  in  the  sample-on-cantilever 
experiment  of  Ref.  25  in  which  4  nm  proton  NMR  imaging  of  a  virus  was  demonstrated. 

The  51  nm  thickness  of  the  magnetically  inactive  layer  in  the  model  is  at  odds  with 
the  <  10  nm  of  oxidation  observed  in  similarly  aged  Co  thin-film  control  samples  using 
depth-resolved  XPS  experiments.  The  apparent  51  nm  thickness  in  the  model  therefore 
likely  arises  from  a  combination  of  oxidation,  leading-edge  roughness,  and  a  protruding 
titanium  adhesion  layer.  This  observation  suggests  that  reducing  the  titanium  overhang 
and  decreasing  the  grain  size  at  the  leading  edge  of  the  magnet  will  reduce  the  thickness 
of  the  apparent  dead  layer,  give  a  larger  tip-field  gradient,  and  improve  the  spin  sensi¬ 
tivity. 

3.2.2  Understanding  non-contact  friction  and  frequency  noise 

When  the  MURI  began,  MRFM  experiments  had  been  carried  out  with  attonewton- 
sensitivity  cantilevers  brought  near  a  surface.  In  the  most  sensitive  of  these  experiments 
it  was  recognized  that  force  sensitivity  was  determined  not  by  the  intrinsic  dissipation 
in  the  cantilever  but  by  the  extra  dissipation  arising  from  cantilever-surface  interactions. 
Rugar  and  co-workers  carried  out  measurements  of  dissipation  over  quartz  and  metal 
surfaces,  showing  that  the  surface  dissipation  arose  from  cantilever  charges  interacting 
with  fluctuating  electric  field  in  the  sample  [17];  over  a  silicon  sample  the  dissipation  was 
shown  to  be  proportional  to  the  sample's  doping  level  [26]. 
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With  partial  support  from  the  MURI  and  primary  support  from  the  NIH  and  NSF, 
Marohn's  team  carried  out  experimental  and  theoretical  work  to  understand  surface  noise 
over  polymer  films  [3-6].  A  first-principles  theory  was  developed  for  electric  field  and 
electric  field  gradient  fluctuations  over  dielectrics,  starting  from  the  stochastic  Maxwell 
equations  of  Lifshitz  and  using  the  sample's  dielectric  function  as  the  only  input  [5].  This 
theory  quantitatively  explained  the  dependence  of  both  dissipation  and  frequency  noise 
as  a  function  of  function  of  polymer  composition,  film  thickness,  and  tip-sample  separa¬ 
tion. 

These  studies  were  applicable  to  a  cantilever  operated  in  the  "hangdown"  geometry 
of  Fig.  1(a).  This  reporting  period,  we  extended  our  theory  to  encompass  a  cantilever 
vibrating  normal  to  the  sample  surface  [7]  —  the  geometry  used  in  some  MRFM  experi¬ 
ments  [27-29]  and  in  essentially  all  atomic  force  microscope  experiments.  This  seemingly 
simple  extension  was  surprisingly  involved  because  now  the  tip  charge  is  time  dependent 
during  a  single  cycle  of  the  cantilever  motion.  We  again  observed  excellent  agreement 
between  theory  and  experiment.  These  findings  open  the  door  to  achieving  quantitive 
dielectric-function  image  contrast  in  AFM  via  measurements  of  dissipation  and  frequency 
noise  versus  height  and  location. 

This  spring  our  theory  collaborators.  Professor  Roger  Loring  and  his  group,  succeeded 
after  many  years  of  effort  in  developing  a  theory  for  frequency  noise  over  a  molecular 
film  containing  mobile  charges.  This  theory,  again  based  on  Maxwell's  equations,  treats 
exactly  the  interactions  of  the  charges  with  each  other  and  with  dielectric  fluctuations  in 
the  host  material.  The  only  inputs  to  the  theory  are  the  material's  dielectric  function,  the 
charge  density,  and  the  charge's  diffusion  constant.  The  theory  makes  the  remarkable 
prediction  that  charge-charge  interactions  suppress  voltage  fluctuations,  and  therefore 
cantilever  frequency  noise,  by  orders  of  magnitude  over  an  organic  transistor  operated  at 
typical  charge  densities  (e.g.,  gate  voltages).  This  prediction  of  suppressed  noise  was  ver¬ 
ified  by  experiment.  A  long  paper  describing  these  theoretical  and  experimental  findings 
has  been  submitted  to  the  Journal  of  Chemical  Physics  [8]. 

It  must  be  kept  in  mind  that  the  sample,  and  not  the  cantilever,  has  set  the  force  sen¬ 
sitivity  in  the  most  sensitive  MRFM  experiments  carried  out  to  date.  Taken  together, 
our  friction  and  frequency  noise  findings  thus  have  important  implications  for  designing 
high-sensitivity  MRFM  experiments. 

Over  dielectric  films  such  as  polymers  and  biopolymers,  our  work  establishes  that 
thermal  dielectric  fluctuations  create  attonewton-level  force  noise  that  can  be  predicted 
from  first  principles.  In  MRFM  experiments  the  sample  is  often  routinely  coated  with  a 
metal  film  to  shield  the  tip  from  the  sample's  electric  field  fluctuations.  In  our  friction  and 
frequency  noise  experiments,  we  find  that  some  polymers  such  as  polystyrene  have  even 
lower  surface  noise  than  gold.  This  finding  suggests  that,  for  some  samples,  gold-coating 
of  samples  is  not  necessarily  the  best  practical  strategy  for  obtaining  low  surface  noise  in 
an  MRFM  experiment.  Our  work  indicates  that,  due  to  the  long  carrier  screening  length, 
charge-charge  interactions  in  molecular  conductors  and  semiconductors  are  effective  at 
suppressing  force  noise  and  frequency  noise.  Thus,  for  "high  noise"  samples,  a  molecular 
conductor  or  semiconductor  coating  may  be  a  preferable  and  lower-noise  coating  than  a 
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metal. 


3.2.3  Spin  detection  protocols 

Before  the  project  began,  we  had  invented  a  force-gradient  method  for  observing  magnetic 
resonance  [30]  —  Cantilever  Enabled  Readout  of  Magnetization  Inversion  Transients  or 
CERMIT  [30;  31].  This  spin-detection  protocol  has  two  unique  capabilities.  The  first  is 
that,  in  principle,  it  requires  only  one  or  a  few  spin  flips  to  observe  signal.  This  capability 
was  recently  exploited  by  Alexson  and  Smith  at  the  U.S.  Army  Research  Laboratory  to 
measure  the  spin-lattice  relaxation  of  69Ga  magnetization  at  various  temperatures  [32]  in  a 
real-time  single-shot  inversion  recovery  experiment.  This  work  was  enabled  by  the  MURI 
because  the  cantilevers  for  the  experiment  were  produced  by  Hickman  and  Marohn  at 
Cornell. 

The  second  advantage  of  the  CERMIT  approach  comes  into  play  when  it  is  used 
with  cantilevers  operating  in  close  proximity  to  a  surface.  High-compliance,  attonewt on- 
sensitivity  cantilevers  must  be  operated  in  the  "hangdown"  geometry  of  Fig.  1(a)  to  avoid 
a  destructive  "snap  in"  to  contact  with  the  surface.  In  this  geometry,  force-based  methods 
for  detecting  magnetic  resonance  can  only  observe  a  lateral  spin  imbalance  in  the  sample, 
such  as  the  imbalance  created  by  magnetization  fluctuations  in  small  ensembles  of  spins. 
The  second  advantage  of  the  CERMIT  approach  is  that  it  enables  the  observation  of  both 
the  mean  spin  polarization  and  polarization  fluctuations. 

With  partial  support  from  the  MURI  and  primary  support  from  the  NIH  and  NSF,  we 
have  continued  to  evaluate  spin-detection  protocols. 

During  the  reporting  period  we  published  a  variant  of  the  CERMIT  approach  that 
extends  the  applicability  of  force-gradient  detection  to  spins  having  spin-lattice  relaxation 
times  Tj  as  short  as  a  single  cantilever  period  of  ~  200  /is  [10].  This  innovation  made  it 
possible  for  us  to  detect,  at  record  sensitivity,  electron  spin  resonance  from  the  stable 
organic  radical  TEMPAMINE  widely  used  in  conventional-ESR  studies  of  proteins  and 
nucleic  acids.  This  radical's  short  spin  relaxation  times  (T\  =  1  ms  and  T2  =  500  ns)  make 
it  impossible  to  observe  using  the  i-OSCAR  force  detection  protocol  employed  by  Rugar 
et  al.  to  observe  single-electron  magnetic  resonance  via  MRFM  [9;  24], 

Our  new  CERMIT  experiment  opens  up  a  route  for  achieving  single-electron  detec¬ 
tion  of  magnetic  resonance  in  organic  radicals  at  reasonably  short  averaging  times.  In 
the  i-OSCAR  experiment,  the  spin  had  a  root-mean-square  magnetization  of  1  /i^;  in  a 
CERMIT  experiment,  one  observes  the  mean  (Curie  law)  spin  magnetization  which  can 
be  as  high  as  0.6 at  2.1  K  and  2.1  T  (e.g.,  at  a  60  GHz  resonance  frequency).  Although  the 
magnitude  of  the  signal  in  a  single-spin  CERMIT  experiment  is  slightly  smaller  than  in  a 
single-spin  i-OSCAR  experiment,  as  discussed  in  Ref.  10  and  shown  below,  the  fact  that 
the  CERMIT  experiment  observes  the  mean  spin  polarization  (which  has  a  well  defined 
sign)  instead  of  polarization  fluctuations  (which  have  random  sign)  makes  the  signal-to- 
noise  in  the  CERMIT  experiment  far  more  amenable  to  improvement  by  signal  averaging. 

A  potential  disadvantage  of  CERMIT  detection  is  that  the  main  source  of  noise  is 
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due  to  surface-induced  frequency  fluctuations,  which  are  in  general  orders  of  magnitude 
larger  than  surface-induced  force  fluctuations.  In  this  reporting  period  we  demonstrated 
using  non-degenerate  parametric  amplification  to  transform  a  slowly  modulated  force- 
gradient  spin  signal  into  a  force  signal  oscillating  on  resonance  with  the  cantilever  [11]. 
This  new  detection  paradigm  has  the  advantages  of  both  force  detection  (less  susceptible 
to  surface  noise)  and  force-gradient  detection  (the  ability  to  observe  both  the  mean  spin 
magnetization  and  magnetization  fluctuations,  potentially  longer  spin-coherence  times 
during  modulation,  and  lower  rf  duty  cycle  and  thererfore  less  heating).  We  demon¬ 
strated  the  approach  by  evading  surface  frequency  noise  and  detector  noise  in  an  MRFM 
measurement  of  ESR  from  TEMPAMINE. 

In  most  prior  theoretical  treatments  of  the  CERMIT  effect,  we  assumed  that  the  can¬ 
tilever  amplitude  was  much  smaller  than  the  diameter  of  the  magnetic  tip  d  and  the 
tip-sample  separation  h  [10;  11;  30;  33].  While  this  approximation  has  been  valid  in  all 
CERMIT  experiments  carried  out  to  date,  optimizing  the  signal-to-noise  ratio  in  a  single¬ 
spin  CERMIT  experiment  requires  considering  the  case  that  the  cantilever  amplitude  is 
comparable  to  d  or  h.  We  considered  this  limit  in  Ref.  33,  but  only  for  a  single  spin 
located  directly  below  the  tip.  This  reporting  period  we  developed  a  formal  description 
of  frequency-shift  MRFM  experiments  valid  at  any  cantilever  amplitude  [34],  Remark¬ 
ably,  we  were  able  to  invent  a  single  unified  approach  to  numerically  calculating  signal 
from  extended  objects  measured  in  both  OSCAR  and  CERMIT  experiments.  The  results 
of  Ref.  34  will  help  us  design  and  simulate  optimal  single-spin,  as  well  as  few-spin,  mea¬ 
surements. 

During  the  course  of  carrying  out  frequency-shift  cantilever  magnetometry  exper¬ 
iments  [21-23]  to  characterize  the  magnetization  of  individual  nanomagnets  on 
attonewton-sensitivity  cantilevers,  we  discovered  a  new  method  for  detecting  nanoscale 
magnetic  fields.  This  discovery  was  published  during  this  reporting  period  [35].  In  these 
experiments  we  applied  a  magnetic  field  oriented  along  the  short  axis  of  the  nanomagnet 
tip  —  in  the  direction  of  the  cantilever  thickness  in  Fig.  1(b).  Since  the  magnetization  is 
oriented  along  the  long  axis  of  the  nanomagnet  at  zero  field,  in  this  experiment  we  are 
observing  a  switching  of  the  nanomagnet's  magnetization  from  its  magnetic  easy  axis  to 
its  magnetic  hard  axis.  This  kind  of  magnetization  reorientation  has  only  been  observed 
once  before  in  a  cantilever  magnetometry  experiment  [36],  and  never  before  with  such  a 
small  magnetic  tip  and  with  such  fine  field  steps.  In  our  experiment  we  observed  multiple 
distinct  sharp  transitions  in  cantilever  frequency,  dissipation,  and  frequency  noise.  These 
features,  which  we  believe  are  due  to  switching  of  individual  magnetic  domains,  are  only 
a  few  gauss  wide.  The  observed  giant  tip-induced  changes  in  cantilever  frequency  and 
dissipation  may  thus  represent  a  new  route  to  detecting  magnetic  fields  (and  ultimately 
magnetic  resonance)  at  nanoscale  spatial  resolution. 

During  this  reporting  period  we  also  investigated  using  dynamic  nuclear  polariza¬ 
tion  (DNP)  [37]  to  align  nuclear  spins  in  an  MRFM  experiment.  Below  we  show  that 
aligning  nuclear  spins  by  DNP  could  dramatically  improve  our  ability  to  detect  and 
image  biomolecules  at  single  proton  sensitivity. 

The  experiment  is  sketched  in  Fig.  2(left).  Experiments  took  place  in  vacuum,  at  a 
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(a)  pulse  sequence 


1 

(c)  prompt 

Figure  2:  Scanned-probe  dynamic  nuclear  polarization  experiment.  Left:  Schematic  drawing  of 
the  experimental  setup.  The  drawing  is  not  to  scale.  The  black  arrow  represents  the  direction  of 
the  longitudinal  polarizing  magnetic  field  Bq.  An  artistically-rendered  spin  indicates  the  sample's 
electron  spin  magnetization,  excited  by  the  Bi.  Right:  (a)  Microwave  pulse  sequence  used  to 
observe  spin-induced  transient  changes  in  the  cantilever's  resonance  frequency;  the  solid  "ON" 
black  bars  indicate  pulses  in  which  the  microwave  frequency  is  in  resonance  with  the  applied  mag¬ 
netic  field.  The  empty  "OFF"  bars  indicate  pulses  in  which  the  microwave  frequency  is  out  of  res¬ 
onance  with  the  applied  field,  (b)  The  observed  cantilever  frequency  shift  as  function  of  time  (red 
circles)  and  fit  to  a  piecewise  sum  of  a  square  wave  and  an  exponential  rise  and  decay  (solid  line). 
Experimental  parameters:  cantilever  amplitude  xc  =  131.7  nm,  tip-sample  separation  h  =  50  nm, 
and  static  field  Bq  =  0.600  T.  (c)  The  prompt,  square-wave  best-fit  component  of  the  observed 
frequency  shift  due  to  saturated  electron  spins  in  the  sample,  (d)  The  prompt  signal  subtracted 
from  the  observed  signal  to  yield  the  persistent  signal  (red  circles).  Also  displayed  is  a  best-fit  to  a 
rising  and  decaying  exponential  (solid  line).  Figure  and  caption  adapted  from  Ref.  38. 


temperature  of  T  =  4.2  K,  and  in  a  background  field  centered  at  Bq  =  0.6  T  corresponding 
to  an  electron-spin  resonance  frequency  of  18  GHz.  A  200  nm  thick  film  of  perdeuterated 
polystyrene  was  doped  with  the  free-radical  TEMPAMINE  at  a  concentration  of  40  mM. 
An  attonewton-sensitivity  cantilever  was  prepared  with  a  4  micron  diameter  nickel  tip 
and  an  optical  fiber  was  aligned  to  the  cantilever  to  watch  its  motion.  The  cantilever  was 
brought  in  close  proximity  to  a  coplanar  waveguide,  which  was  shorted  at  the  end  to 
generate  a  large  transverse  microwave  magnetic  field  Bi  (the  broad  arrow  in  Fig.  2(left)). 
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The  cantilever  was  parked  above  the  middle  of  one  of  the  shorted  segments.  Changes  in 
the  sample  magnetization  were  registered  as  a  shift  in  the  mechanical  frequency  of  the 
cantilever  [10]. 

The  time  sequence  of  applied  microwaves  and  the  resulting  cantilever  frequency  shift 
are  plotted  in  Figure  2(right).  Microwaves  were  applied  continuously  and  their  frequency 
was  switched  essentially  instantaneously  from  off-resonance  to  on-resonance.  Because 
the  applied  power  was  constant,  shifts  in  cantilever  frequency  induced  by  switching- 
related  thermal  transients  were  negligibly  small.  Upon  shifting  the  microwave  irradia¬ 
tion  on  resonance,  the  cantilever  frequency  exhibited  a  prompt,  essentially  instantaneous 
change  due  to  saturation  of  the  unpaired  electron  spins  in  the  TEMPAMINE  radical. 
Superimposed  on  this  instantaneous  frequency  shift  was  a  slower,  persistent  frequency 
shift.  The  risetime  of  the  persistent  shift  was  much  longer  than  any  thermal  transient.  As 
the  applied  field  B0  was  scanned,  the  magnitude  and  sign  of  the  persistent  shift  tracked 
the  magnitude  and  sign  of  the  prompt  shift  (data  not  shown).  Because  of  its  timescale, 
dependence  on  applied  field  B0/  magnitude,  and  exponential  kinetics,  we  attribute  the 
persistent  shift  to  buildup  of  nuclear  magnetization.  A  paper  describing  these  findings 
has  been  submitted  [38]. 

A  persistent  shift  was  not  observed  when  the  applied  field  was  set  in  resonance  with 
local  spins  directly  below  the  tip;  buildup  was  only  observed  when  the  static  field  was  set 
in  resonance  with  bulk  spins  far  away  from  the  tip.  This  finding  indicates  that  buildup 
of  dynamic  nuclear  polarization  of  deuterium  may  be  suppressed  by  the  large  magnetic 
field  gradient  present  directly  below  the  tip.  Flarnessing  dynamic  nuclear  polarization 
in  an  MRFM  experiment  may  therefore  require  minimizing  the  magnetic  field  gradient 
experienced  by  the  spins  during  DNP.  We  propose  to  achieve  this  in  future  experiments 
by  shuttling  the  cantilever  to  the  side  during  the  DNP  period. 

3.2.4  Image  Reconstruction  Algorithms 

We  have  focused  on  building  models  and  algorithms  for  3D  MRFM  image  reconstruction 
when  the  point  spread  function  is  uncertain  or  only  partially  known.  This  problem  is 
significant  and  our  models  and  approaches  have  been  strongly  influenced  by  discussions 
with  Dan  Rugar. 

During  this  period  we  have  made  significant  progress.  We  published  journal  papers 
in  IEEE  Trans,  on  Image  Processing  in  Sept.  2009  [39]  and  Sept.  2012  [40],  and  one  other 
paper  is  in  review  at  the  same  journal  [41].  The  2009  journal  paper  introduced  a  hierar¬ 
chical  Bayesian  method  for  3D  image  reconstruction  when  the  point  spread  function  is 
known  exactly.  The  method  provides  a  comprehensive  approach  to  MRFM  reconstruc¬ 
tion  that  incorporates  physics  models  of  the  image  formation  process  and,  in  addition 
to  the  image  reconstruction,  provides  confidence  intervals  on  the  voxel  intensities  of  the 
reconstructed  image.  Such  confidence  intervals  are  essential  for  doing  hypothesis  testing 
on  reconstructed  images,  e.g.  detecting  an  edge  or  shape  of  an  object  in  the  image,  and 
were  not  available  using  standard  Landweber  methods  used  by  Rugar.  The  method  was 
applied  to  image  reconstruction  of  the  tobacco  virus  data  collected  by  Rugar's  group  at 
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IBM. 

The  two  other  papers  significantly  extend  the  work  described  in  the  2009  IEEE  paper 
to  the  case  where  the  point  spread  function  is  uncertain.  The  first  paper  reports  on  a 
MCMC  approach  to  estimating  simultaneously  the  point  spread  function  and  the  image 
intensities,  solving  the  so-called  semi-blind  reconstruction  problem.  The  second  paper 
introduces  a  variational  Bayes  approach  that  implements  the  semi-blind  reconstruction 
approach  of  the  first  paper.  This  variational  Bayes  approach  has  significant  computation 
advantages.  Both  semi-blind  approaches  come  with  confidence  intervals  on  image  inten¬ 
sities  and  on  point  spread  function  parameters. 

In  summary,  we  have  developed  a  new  and  flexible  class  of  MRFM  image  recon¬ 
struction  methods.  The  advantages  of  these  methods  are  the  following.  Unlike  previous 
methods, 

•  the  approach  directly  imposes  natural  sparsity  models  on  the  image, 

•  there  are  no  sensitive  tuning  parameters  to  select, 

•  the  reconstruction  incorporates  uncertainty  on  point  spread  response, 

•  the  reconstruction  algorithm  outputs  a  full  probability  distribution  of  each  pixel 
intensity, 

•  the  performance  of  the  algorithm  can  be  monitored  on  the  basis  of  this  distribution, 
and 

•  the  algorithm  converges  to  the  global  maximum  of  the  Bayesian  risk  function  asso¬ 
ciated  with  mean  square  reconstruction  error. 


3.2.5  Spin  Simulation  Algorithms 

Simulations  play  a  central  role  in  modern  science  and  engineering,  but  it  remains  a  chal¬ 
lenge  to  simulate  the  large-scale,  high-temperature,  and  nonequilibrium  quantum  sys¬ 
tems  that  occur  in  applications  such  as  magnetic  resonance  imaging,  condensed  matter, 
chemistry,  and  nanotechnology.  Geometric  descriptions  of  quantum  trajectory  unraveling 
lead  to  algorithms  that  can  speed  the  computation  of  quantum  simulations  for  broad 
classes  of  physical  systems1. 

Simulation  speed  is  enhanced  whenever  the  rank-n  matrix-vector  products  that  gener- 
ically  appear  in  simulation  equations  are  structured  so  as  to  be  computed  in  O(n)  oper¬ 
ations,  rather  than  the  0(n 2)  generally  required.  In  recent  years  the  resulting  n-fold 
speedup  has  transformationally  augmented  the  feasible  scale  of  classical  simulations. 

During  this  reporting  period  we  have  been  concerned  with  the  question:  Can  a  sim¬ 
ilar  n-fold  speedup  be  achieved  in  the  simulation  of  quantum  spin  systems?  We  find  that 

1The  text  in  this  section  is  drawn  nearly  verbatim  from  Ref.  42 
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the  answer  is  "yes."  A  long-established  principle  of  of  quantum  system  analysis  simu¬ 
lation  is  that  "a  very  large  number  of  different  physical  viewpoints  and  widely  different 
mathematical  formulations  that  are  all  equivalent  to  one  another."  We  have  exploited  this 
freedom  to  develop  a  quantum  simulation  framework  that  is  closely  similar  to  the  sym- 
plectic  mathematical  framework  of  classical  molecular  dynamical  simulation,  and  that 
simultaneously  respects  the  orthodox  principles  of  quantum  mechanics,  as  these  princi¬ 
ples  appear  when  translated  onto  symplectic  state-spaces. 

The  central  idea  is  to  implement  the  orthodox  equations  of  quantum  simulation  in 
such  a  way  that  simulated  trajectories  do  not  fill  the  entire  Hilbert  space,  but  rather  are 
dynamically  compressed  onto  a  curved  state-space  of  lower  dimension,  such  that  the 
curved  state-space  geometry  describes  the  trajectory  implicitly  rather  than  explicitly.  This 
is  a  familiar  strategy  in  simulating  classical  systems  by  projective  model  order  reduction 
(MOR),  but  it  is  less  widely  appreciated  that  this  strategy  is  similarly  effective  in  quantum 
simulations. 

As  an  illustrative  case  we  have  simulated  a  generic  dynamic  nuclear  polarization 
(DNP)  process,  in  which  spin  polarization  is  transferred  from  a  strongly  polarized  elec¬ 
tron  spin  to  an  unpolarized  nuclear  spin. 

It  is  characteristic  of  fields  like  biology  and  materials  science  that  large  ensembles  of 
spins  (or  atoms,  or  electrons)  are  not  the  systems  of  interest,  but  rather  investigations 
focus  on  nanoscale  systems  having  a  countable  number  of  particles.  In  such  cases  the 
mean  values  of  quantities  like  spin  polarization  can  still  be  determined  by  taking  time 
averages,  but  fluctuations  in  state-space  processes  and  measurement  processes  are  also 
observed.  Therefore  our  framework  encompasses  the  fluctuations  in  both  the  spin  polar¬ 
ization  and  in  a  measurement  process  that  continuously  monitors  that  polarization,  using 
algorithms  that  scale  feasibly  to  the  simulation  of  systems  of  hundreds,  or  potentially 
even  thousands,  of  dynamical  quantum  spins. 

3.3  Single  proton  detection  outlook 

In  summary,  our  MURI  team  made  tremendous  advances  in 

•  developing  nanomagnet-tipped  attonewton-sensitivity  cantilevers, 

•  microscopically  and  spectroscopically  quantifying  the  chemical  composition  of  pro¬ 
cessed  magnetic  films  and  individual  nanomagnets, 

•  understanding  the  force  and  force-gradient  fluctuations  experienced  by  an  ultrasen¬ 
sitive  cantilever  near  a  surface, 

•  inventing  spin-readout  protocols  that  greatly  extend  the  applicability  of  magnetic 
resonance  force  microscopy, 

•  inventing  spin-readout  protocols  that  evade  surface  frequency  noise  and  detector 
noise. 
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•  developing  a  new  and  flexible  class  of  MRFM  image  reconstruction  methods,  and 

•  inventing  and  reducing  to  practice  new  approaches  for  simulating  many  coupled 
spins. 

In  addition,  we  have  presented  preliminary  data  suggesting,  for  the  first  time  in  an  MRFM 
experiment. 


•  alignment  of  nuclear  spins  by  transfer  of  polarization  from  electron  spins. 


Although  we  did  not  reach  our  goal  of  detecting  single  proton  magnetic  resonance,  our 
MURI  team's  work  has  given  us  a  very  clear  understanding  of  the  remaining  technical 
hurdles  that  need  to  be  overcome  to  achieve  single  proton  sensitivity.  Moreover,  our 
work  suggests  a  clear  technical  path  to  this  goal,  which  we  describe  in  this  section. 

We  begin  by  assessing  the  signal-to-noise  of  the  Cornell/IBM  experiment  described 
above,  and  show  that  single  proton  sensitivity  can  be  reached  from  this  starting  point  via 
a  combination  of  (1)  challenging  but  achievable  reductions  in  the  thickness  of  the  tip's 
magnetically  inactive  layer  and  (2)  alignment  of  nuclear  spins  via  transfer  of  polarization 
from  electron  spins. 

We  begin  by  calculating  the  signal-to-noise  ratio  for  observing  an  individual  proton 
using  the  Cornell/IBM  experiment  of  Fig.  1.  The  signal-to-noise  ratio  for  detecting  mag¬ 
netic  resonance  from  a  stochastically  polarized  single  proton  in  a  force-based  magnetic 
resonance  experiment  is 

SNRp  =  yjT^  rm  (1) 

where  yp  =  1.41  x  10  26  N  m  T  1  is  the  proton  magnetic  moment,  Gzx  =  f)Bfp/ dx  is  the 
lateral  tip-field  gradient,  Tavg  is  the  signal-averaging  time,  r„,  is  the  spin  dephasing  time 
during  modulation,  and  SF  is  the  spectral  density  of  force  fluctuations  experienced  by  the 
cantilever. 


To  estimate  the  gradient  achievable  with  improved  tips,  is  it  convenient  to  model  the 
tip  as  a  uniformly  magnetized  sphere.  The  simplest  tip  model  is  shown  in  Fig.  3(a).  To 
account  for  the  magnetically  inactive  layer  observed  in  all  of  our  tips  to  date,  we  will  use 
the  tip  model  shown  in  Fig.  3(b).  For  comparison,  the  tip  model  used  to  analyze  the  data 
in  Fig.  1  is  shown  in  Fig.  3(c). 


To  magnetize  the  tip,  a  polarizing  magnetic  field  is  applied  along  the  2  axis  in  Fig.  3. 
We  assume,  for  simplicity,  that  the  magnitude  of  the  polarizing  field  is  larger  than  the  tip 
field;  in  this  case  we  only  need  to  consider  the  2  component  of  the  tip  field  in  the  following 
calculations.  The  2  component  of  the  magnetic  field  outside  of  the  uniformly  magnetized 
sphere  is 


(2) 


with  r  =  y/x 2  +y2  +  z2  the  distance  from  the  center  of  the  sphere  to  the  point  of  interest 
r  =  (. x ,  y,  2).  Here  Bs  =  2/i0Msat/3  is  the  magnetic  field  at  the  north  pole  of  the  sphere. 
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Figure  3:  Tip  models  used  to  estimate  the  single-proton  signal,  (a)  Spherical  tip  of  radius  a.  The 
physical  tip-sample  separation  is  h„hys-  A  proton  is  located  a  vertical  distance  2  from  the  center  of 
the  tip.  (b)  Spherical  tip  including  a  (shaded)  magnetically  inactive  layer  of  thickness  /2 inactive-  The 
effective  magnetic  spacing,  the  physical  spacing  plus  the  thickness  of  the  inactive  layer,  is  indi¬ 
cated  as  /imag.  (c)  For  comparison,  the  tip  model  is  shown  that  was  used  to  numerically  simulate 
the  spin  variances  in  Fig.  1(c). 


with  Msat  the  saturation  magnetization  and  Ho  =  47r  x  10  'TmA  1  the  permeability  of 
free  space.  For  cobalt,  /U0Msat  =  1.8  T  and  Bs  =  1.2  T. 

We  need  to  consider  two  field  gradient  components.  The  first  is  the  vertical  tip-field 
gradient  Gzz  =  dB^p/dz.  Directly  below  the  tip,  r  =  (0,  0,  z)  and  GA(0, 0,  z)  =  —  3Bsa3/z4, 
with  z  =  a  +  hmag  (Fig.  3(b)).  At  a  fixed  magnetic  separation  hma&  the  vertical  gradient 
is  maximized  with  a  magnetic  tip  radius  of  aopt  =  3  hmag.  The  vertical  gradient  from  this 
optimized  magnetic  tip  is  at  a  maximum  directly  below  the  tip  and  equal  to 

G°pt  =  -0.316  — .  (3) 

hmag 

The  second  component  to  consider  is  the  lateral  tip-field  gradient  Gzx  =  0B)'P /dx.  This 
gradient  component  determines  the  force  on  the  sample  spins;  see  Eq.  1  above  and  Eq.  5 
below.  For  spins  in  an  x ,  y  plane  located  a  distance  (0,  0,  z)  below  the  tip  center,  the  lateral 
gradient  component  is  maximized  off  to  the  side,  at  locations  x  =  ±0.3892  and  y  =  0.  Sub¬ 
stituting  a  =  aopt  =  3  hmag  gives  the  tip-optimized  lateral  field  gradient  at  these  locations 
as 

G2*  =  ±  0.144  Tk  (4) 

^mag 

By  comparing  Eq.  3  to  Eq.  4,  we  see  that  for  a  spherical  tip  the  tip-optimized  lateral  gra¬ 
dient  off  to  the  side  is  a  factor  of  0.457  smaller  than  the  tip-optimized  vertical  gradient 
directly  below  the  tip. 

The  spherical-tip  model  is  a  reasonably  good  representation  of  the  real  cuboid  tip. 
From  numerical  models  of  the  Fig.  1  signal  we  estimated  the  tip  gradients  in  the  Cor- 
nell/IBM  experiment  to  be  Gzz  =  4.4  to  5.4  MT  m”1  and  Gzx  =  2.7  to  8.3  MT  m_1  at  a 
tip-sample  separation  of  hp hys  =  13.1  nm.  The  lower  bound  is  for  a  model  assuming  a 
cuboid  tip  have  a  leading-edge  damage  layer  of  /inactive  =  51  nm  while  the  upper  bound  is 
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for  a  model  in  which  the  tip  has  uniformly  lower  magnetization  but  no  inactive  layer.  We 
model  the  tip  as  a  sphere  using  /inactive  =  51  nm  (Fig.  3(b))  and  adjust  the  tip  radius  to  give 
Gzz  =  5.4  MT  m-1.  The  resulting  effective  radius  is  a  =  120  nm,  very  close  to  the  133  nm 
geometric  mean  of  the  tip  thickness  (79  nm)  and  width  (225  nm).  The  resulting  Gzx  is 
2.5  MT  m_1,  very  close  to  the  lateral  tip-field  gradient  reported  in  Ref.  20.  This  agreement 
establishes  that  we  can  use  the  spherical  tip  model  to  estimate  the  improvement  in  lateral 
gradient  obtainable  by  decreasing  the  thickness  of  the  inactive  layer. 

We  use  Eq.  1  to  calculate  the  signal-to-noise  for  detecting  nuclear  magnetic  reso¬ 
nance  from  a  stochastically  polarized  single  proton  using  the  magnet-on-cantilever  Cor- 
nell/IBM  experiment  of  Fig.  1  [20].  For  consistency,  we  will  use  the  gradient  computed 
from  the  spherical-tip  model.  The  result  is 

Case  IA.  I  Gzx  =  2.5  MT  m~4 
Cornell  cantilever  now  SF  =  (40  aN)2  FIz”1 
Stochastic  proton  polarization  rm  =  0.1  s 
SNRp  =  1.2  x  10“4  Tavg  =  3  d 

For  comparison,  for  the  high-frequency  cantilever  in  the  sample-on-cantilever  experiment 
of  Nichol,  Budakian,  and  coworkers  [43], 


Case  IB.  I  Gzx  =  0.12  MT  m-1 
UIUC  cantilever  now  SF  =  (1.9  aN)2  FIz-1 
Stochastic  proton  polarization  rm=  1.0  s 
SNRp  =  3.8  x  10~4  Tavg  =  3  d 


which  is  slightly  better  but  SNR  is  still  <C  1.  The  silicon  nanowire  MRFM  experiment 
of  Ref.  43  had  impressive  force  sensitivity,  but  their  detection  method  requires  a  time- 
dependent  gradient  produced  by  running  current  through  a  wire  and  requires  the  sample 
to  be  affixed  to  the  cantilever. 

In  the  magnet-on-cantilever  Cornell /IBM  experiment  we  believe  it  is  feasible  to  both 
increase  the  gradient  and  decrease  the  surface  force  noise.  A  tip  similar  to  that  shown 
Fig.  1(a)  [1]  experienced  a  spectral  density  of  surface  noise  as  good  as  10  aN  FIz”1,72  down 
to  tip-sample  separations  as  small  as  h  =  3  nm  when  operated  over  a  gold-coated  sample 
with  the  tip  voltage  adjusted  to  minimize  the  noise.  Based  on  the  known  thermal  oxi¬ 
dation  layer  thickness  in  cobalt,  corroborated  by  our  own  XPS  depth  profiling  studies  of 
film  oxidation  and  capping,  we  believe  that  reducing  the  thickness  of  the  inactive  layer 
to  ^active  =  5  nm  will  ultimately  be  possible.  At  a  tip-sample  separation  of  hp hys  =  10  nm 
and  /rmag  =  15  nm,  the  optimal  tip  radius  is  aopt  =  45  nm,  and  the  lateral  gradient  is 
Gzx  =  11.5  MT  m_1.  Under  these  conditions, 

Case  IC.  I  Gzx  =  11.5  MT  m_1 
Cornell  cantilever  improved  SF  =  (10  aN)2  FIz-1 
Stochastic  proton  polarization  rm  =  0.5  s 
SNRp  =  0.095  |  Tavg  =  3  d 
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The  signal  to  noise  for  the  stochastically-polarized  proton  experiment,  although  signifi¬ 
cantly  improved,  is  still  well  below  1.  To  obtain  a  signal-to-noise  of  SNRp  =  2.4  would 
require  SF  =  (2  aN)2  Hz-1,  which  has  never  been  achieved  with  an  audio-frequency  can¬ 
tilever  near  a  surface.  It  is  likewise  difficult  for  us  to  see  how  to  improve  the  gradient  in 
the  Ref.  43  experiment  enough  to  obtain  SNRp  >  1  in  reasonable  averaging  times. 

Now  let  us  examine  how  the  signal-to-noise  ratio  improves  by  aligning  the  proton 
magnetization  using  DNP  and  detecting  a  coherent  signal  instead  of  an  incoherent  signal. 
The  signal-to-noise  ratio  for  detecting  magnetic  resonance  from  a  coherent,  well-aligned 
single  proton  in  a  force  experiment  is 

SNRa  =  ^  (5) 

with  fP  <  1  the  spin  polarization  of  the  individual  proton  and  with  G,  fjp,  Tavg,  and  SF 
defined  as  in  Eq.  1.  Here  we  assume,  for  simplicity,  that  the  time  to  polarize  the  proton  is 
shorter  than  the  spin  coherence  time  during  detection;  that  is,  we  neglect  any  "deadtime" 
penalty  associated  with  polarizing  the  proton  magnetization.  We  expect  this  penalty  to  be 
small  since  both  the  DNP  time  and  the  spin  coherence  time  during  detection  are  thought 
to  approach  the  nuclear  spin's  Tx  with  strong  rf  and  large  frequency  modulations. 

In  a  DNP  experiment  the  proton  can  in  principle  achieve  a  polarization  fp  equal  to  the 
electron  spin  polarization  fe,  given  by 


f e  =  tanh  (  Te  °)  (6) 

where  h  =  1.054  x  10~34  J  s_1  is  Planck's  constant,  ye  =  1.76  x  1011  rad  s_1  T”1  is  the 
gyromagnetic  ration  of  the  electron,  /rB  =  1.38  x  10“  23  J  K“ 1  is  Boltzmann's  constant,  B0 
is  the  applied  field  in  tesla,  and  T  is  the  temperature  in  kevin.  In  Fig.  4  we  plot  the 
electron  spin  polarization  achievable  at  T  =  0.3  K  (a  pumped  3He  refrigerator),  T  = 
2.1  K  (pumped  4He),  and  T  =  4.2  K  (liquid  helium)  as  functions  of  both  magnetic  field 
and  the  corresponding  electron  spin  resonance  frequency.  At  T  =  2.1  K  and  B0  =  2.1  T 
(e.g.,  60  GHz),  the  calculated  electron  spin  polarization  is  pe  =  0.60.  Assuming  perfectly 
efficient  polarization  transfer  from  electrons  to  protons,  fp  =  0.60  also.  Using  the  Cornell 
tip  of  Fig.  1,  the  signal-to-noise  for  detecting  a  coherent  single  proton  would  be 

Case  IIA.  I  Gzx  =  2.5  MT  m~4 
Cornell  cantilever  now  SF  =  (40  aN)2  Hz-1 
Aligned  proton  polarization  fp  =  0.60 
SNRa  =  0.27  |  Tavg  =  3  d 

which  is  a  three  order-of-magnitude  improvement  over  the  incoherent-spin  force  exper¬ 
iment  of  Case  IA.  Improving  the  force  noise  to  the  Hickman  et  al.  value  of  10  aN  Hz  1/2 
improves  the  signal-to-noise  to  SNRA  =  1.0  —  quite  exciting.  Interestingly,  the  estimated 
signal-to-noise  for  an  analogous  coherent  single-proton  experiment  carried  out  the  UIUC 
cantilever  is  similar: 
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electron  spin  resonance  frequency  [GHz] 


magnetic  field  [T] 


Figure  4:  Electron  spin  polarization  fe,  Eq.  6,  versus  applied  magnetic  field  Bq  in  tesla  (lower  axis). 
The  corresponding  electron  spin  resonance  frequency  fe  =  'yeBo/2ir  in  GHz  is  shown  on  the  upper 
axis.  The  electron  spin  polarization  is  shown  at  a  temperature  T  of  4.2  K  (red),  2.1  K  (black),  and 
0.3  K  (blue).  The  solid  black  horizontal  line  indicates  an  electron  spin  polarization  of  pe  =  0.60. 


Case  IIB.  I  Gzx  =  0.12  MT  m-1 
UIUC  cantilever  now  SF  =  (1.9aN)2Hz_1 
Aligned  proton  polarization  fp  =  0.60 
SNRA  =  0.26  |  Tavg  =  3  d 


As  long  as  fp  >  0.003  =  0.03%,  the  signal  to  noise  for  the  coherent  experiment  exceeds  that 
for  the  incoherent  experiment.  Assuming  well-polarized  protons  in  combination  with  the 
reduced  noise  and  improved  tip  used  in  Case  IC  above,  we  calculate 


Case  IIC. 

Cornell  cantilever  improved 
Aligned  proton  polarization 
SNRa  =  4.95 


Gzx  =  11.5  MT  m  1 
SF  =  (10  aN)2  Hz’1 
fP  =  0.60 

-Gvg  =  3  d 


Relaxing  the  signal  averaging  time  to  1  day  gives  SNRa  =  2.85  which  is  enough 
(assuming  Gaussian  noise)  to  identify  more  than  99%  of  protons  in  an  individual 
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molecule. 

We  believe  the  parameters  used  in  the  Case  IIC  calculation  are  challenging  but  achiev¬ 
able.  Already  we  have  evidence  of  <  5nm  of  surface  damage  in  our  best  processed  cobalt 
films,  comparable  to  the  thermal  oxidation  thickness.  DNP  at  30%  efficiency  is  routine 
and  DNP  on  the  seconds  timescale  has  been  achieved  in  conventional  magnetic  resonance 
experiments  [37],  though  not  a  temperatures  as  low  as  2.1  K. 
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